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^ '■ Abstract 

CS I The evidently supersymmetric four-dimensional Wess-Zumino model with quenched 

j~^ ■ disorder is considered at the one-loop level. The infrared fixed points of a beta- 

ON ■ function form the moduli space Ai = RP^ where two types of phases were found: 

with and without replica symmetry. While the former phase possesses only a 

trivial fixed point, this point become unstable in the latter phase which may be 

^' interpreted as a spin glass phase. 

> 
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The Devil Is Not So Black As He Is Painted. 



1 Introduction 

There is a great deal of field theory models describing a system in quenched ran- 
dom fields or coupling constants ([0], p|, 0, etc.). In solid state physics such models 
naturally arise the corresponding pure systems whenever impurities are introduced. It 
is interesting to extend randomness to other well-studied field theories, just as, for ex- 
ample in [P], disorder was implemented into minimal conformal models. As shown in [Q 
and subsequent papers stochastic equations as well as the field theories in presence of 
random external sources often prove to possess some hidden supersymmetry. Kurchan 
indorsed this result for spin glass dynamics. Because supersymmetry handle pertur- 
bative corrections, such random theories are especially interesting. This is what we will 
do in this paper. 

On the other hand in field theories with apparent space-time supersymmetry super- 
potential is holomorphic function not only of fields but also of coupling constants [0. 
Therefore couplings and fields enter potential on equal footing, so that it seems very nat- 
ural to introduce random (gaussian) distribution of some couplings in the Lagrangian. 
But the power of supersymmetry is so strong that superpotential gets no quantum cor- 
rections 0,0, i-e. provided that the coupling has no dynamical D-terms integrating it 
over solves the problem. 

In Section 2 we formulate four- dimensional supersymmetric Wess-Zumino theory in 
random field. In the context of replica method infrared fixed points of one-loop /3- 
functions are found in Section 3. Analysis of these fixed points suggests two phases 
on the moduli space Ai = RP^. Numerical evaluation of the most general expressions 
is eventuated in the phase diagram which is illustrated by two simple examples of the 
fourth section. Section 5 is devoted to discussions and conclusions. 

2 Wess-Zumino model perturbed by randomness 

From above arguments it follows that the SUSY analog of a theory with disorder 
must contain dynamical terms for the random field. In the present paper we consider 
a four-dimensional Wess-Zumino model that is the supersymmetric counterpart of f^- 
model (the both theories are defined in the same critical dimension and the scalar poten- 
tial after integrating auxiliary field in the former model is actually y^^). Since, according 
to [0, Wess-Zumino theory is defined only as a low-energy field theory, we will study 
the Wilsonian effective action by integrating fast modes with momentum A' < P < A. 
Thereby, let us define a chiral superfield ^ = (p + Oip + 9^F and a random superfield H. 
In this notations the original action is^: 

S= I d^xd^e(fe{g^+^ - ^+H - //+$ + -H+H) + 
J u 



"'^For the sake of simplicity the mass term.s are om.itted. 



+^ / d^x(fe{>!^^H^ + X'^^^H + X'^^"^ + A;//^) + /i.e. (1) 

This action admits the following treatment. It may be obtained (for the certain set 
of parameters) from the usual Wess-Zumino action by the replacement $ ^ $ + i^, as 
one usually does in summation over local extremes |^. 

One of the most powerful method to deal with random fields is the replica trick [|1|, 
which we will use here to solve this "toy" model. It reduces to introducing n copies 
(replicas) of our system, integrating H field out, then solving n-replica problem and 
taking n = at the end of calculations. After replication the action ([^) takes the form : 



•' a=l "^ 

1 /■ " 

+ / d'xd%Y.(x:,^aH^ + A^$2// + a;,*:^) + x:,h'] + h.c. 



(2) 



As will be shown later the model depends only on the relative values of lambdas, so 
that one can put them small enough to determine H field from the saddle point equation 
on D-term only: 

n n 

E = uY,^a and H+ = uY,^l (3) 

a=l a=l 

Substituting it back into (H) yields: 

S = E:t=i I d'xd'ed'egab^t^t + (,. 

H, I d'xd^9{El,,,^, Ai$„$,$, + Ea%=i A2$I$6 + ELi A3$D + h.c. ^^^ 

where gaa = 9 + 3m, gaj^b = 3m and three types of vertexes Ai = X[u^ + X'^^u^, A2 = AgW, 
A3 = A3 band differently replica indices. It is the action (^ that we are going to study. 

3 Fixed points of /^-functions 

Renormalisation group (RG) equations for gab easily follow from the one-loop dia- 
gram for the pure Wess-Zumino theory P] : 

dn 1 " 

— = ^{SAg^^fc + 2A2 J2 iidac + 9bc)9cd + 9ac9bd] + 



d\nK 28877 , -, 

c,a=l 



+3A2A3 



Y.^9lc + 9lc) + 2^ab Y.(9ac + gbc 
.c=l c=l 



^>^lh ^ {.9ac9ad + gbcgbd)} (5) 

c4=l 



Taking into account the possible replica symmetry breaking we take the Parisi ansatz 
for gab [0: off-diagonal part of gab is parametrized by internal function g{x) defined on a 
unite interval x G [0, 1] and diagonal part is gaa = 9. Replica-symmetric case is obtained 
by putting g{x) = g = constant. Algebra of Parisi matrices a = {d,a{x)) is defined by 




g« 



+-^C:> 



I > 



Figure 1: Surviving (in the n ^ limit) A2-contributions. 

the multiphcation rule |]1|: 

c = ab : c= ah — I dxa{x)b{x) 

Jo 

c{x) = b{x)[a — / dxa{y)] + a{x)\b — / dx6(y)] — (6) 

JO Jo 

dy{a{x) - a{y)){b{x) - h{y)) 

By means of this rule we get sums over replica indices that appear in (^ in the n -^ 
limit: 

n n n 

Y.9ac^g-g Y. dacQcd ^ {g- gf Y. die ^ f - g^ (7) 

6=1 c,d=l 6=1 

where 

g = dxg{x) and g^ = dxg'^{x) (8) 

JO JO 

A question arises, as usual in spin glass theory, to find infrared (IR) fixed points^ of 
(^ which determine the dynamics of the system: 

\\lf + (A^ + 3AiA3)(^ - -gf + A2A3[2^(^ -g)+g^- ?] = (9) 

\\lg\x) + (A^ + 3AiA3)(^ - gf + \2H2g{x){g --g)+f- ?] = (10) 

For example, the Al-term is produced by the two non-vanishing (with the number of 
replicas) diagrams shown on Fig.l. 

These equations have two remarkable properties: they are homogeneous in A and g, 
i.e. depend only on the squares of the both. Such dependence on A tells us that zeroes 
of beta- functions (pD-(p!OD do not depend on the values of the couphngs themselves, but 
only on their mutual ratios, so that the moduli space of the theory is RP^ instead of 
E? = {Ai, A2, A3}. Therefore, without loss of generality, we may put couplings very small 
keeping their ratios fixed. In this limit the results that we are going to obtain are exact. 
Moreover, in what follows we will assume A3 7^ 0, so that we can choose it to be A3 = 1 
and denote A2 = A and Ai = /i (affined map). | The special case of A3 = will be 
studied in the first example of Section 4. 

^The points where /3-functions vanish. 

•^If A3 7^ 1 then the right parameters are ^ — j^ and fi — j^. 



Quadratic dependence on g in (|TU|) means that for each set of general characteristics, 
such as g, g'^ and g, there are only two possible values gi^2 (if any) which the function g{x) 
can take in a IR-fixed point. Moreover, the same must be true for g because formally 
it also satisfies similar equation (||). We are free to chose g = gi, for instance. Let us 
denote the measure of points on a unite interval of x where g{x) = gi as 1 — Xq and the 
measure of points where g{x) = g2 as Xq. For example, it may be a stepwise distribution: 

Thus we have two equations (^-(0) o^ three quantities gi^2 and xq with g and g"^ 
depending on them. If gi and g2 are not simultaneously equal to zeroQ then, actually, 
we have only two unknowns: Xq and the ratio p = —. In this notations (0)-(0) niay be 
rewritten as: 

A^ + 2/i)x2(l - pf + |xoA [2(1 -p) + {l- p')] = 

A2 + 2/i)x2(l - p)2 + |xoA [2p{l - p) + (1 - p^)] = ^'^> 

determining both p and Xq, and, consequently, the phase of the system. 

Curiously enough, for a given solution p and Xq we get the whole set of RG-fixed 
points {g,g{x)} differing in arbitrary factor. Of course, this degeneracy will be removed 
by higher loop corrections, so that particular value of the fixed point will be determined 
by the full perturbative expansion. At the one-loop approximation, the explicit data 
{g,g{x)) in the fixed point may be determined by the initial conditions g and u. 

If for some set of couplings there is no solution to (p^ ) except the trivial one g = 
g{x) = 0, we will refer to this point on the phase space {A, fi} & Ai = RP^ as a replica- 
symmetric point and will denote the corresponding phase "RS". Otherwise, replica 
symmetry is broken with xq being the solution of (|1^), and the corresponding phase 
" RSB" looks like a spin glass system. 

Since (|1^) must be solved by the same p, equating the solutions to each equation we 
get the relation between xq and {A, //} G Ai. Instead of writing the resulting complicated 
formulae (partly because it can not be resolved relatively Xq), we display it for xq = 1: 




A2 + 3^ + A ± y'f A2 - I;, + |A _ A2 + 3/i - A ± ^fA^ - |/i - |A 
A2 + 3/i - A ~ I + A2 + 3/i - 3A 



(13) 



where sings in the both parts are taken independently. Replacing X —>■ xqX and /i -^ x^fi, 
we get the equation ([T3|) for arbitrary xq. This expression describes (part of) a curve in Ai 
that separates RS and RSB phases as shown on Fig. 2. The shaded region indicate replica- 
symmetric phase and the unshaded region corresponds to replica symmetry breaking 
where there is a non-trivial solution to ([T3|) , and the trivial point g = g{x) = becomes 
unstable, as will be descanted in the second example of the next section. 

^Otherwise we get a trivial replica-symmetric fixed point. 




Figure 2: The Phase Diagram (not drawn in scale). 



4 Two simple examples 

• A3 = 

In this case the beta-functions (^) become 



dg 



din A 

dg{x) 



487r2 
1 



d In A 4871^ 



m-gr 

\li~9-9f 



These equations may be easily integrated with the result: 



A 



A 



9k{.x)= ^o,A' {.x) + — ^ A2 In — 



487r2 



A' 



where a constant A= {g — gY is determined by initial conditions and remains un- 
changed during renormalisation group flow. Since for any A2 the only fixed point 
is ^ = g{x) = 0, this phase is always replica-symmetric and is not as interesting as 
others. 



Aa = *-* A = 

Equations (|9|)-([10|) take the form: 



g\x)+2ii{~g--gf = Q 



(14) 



for which g'l 2 = ±(7 for some (7 7^ in the SG phase. In parametrization (|Tl|) 

9 ={.91- 92)xo = 2gxo and f = {gf - gl)xo = 2g^xo (15) 

Substitution it into ( [I^ ) yields a nontrivial solution: 

- S/iXg = 1 or Xq = -^=== (16) 

which exists only for // < — |. It is the range of /i where the RSB phase can be 
found. Let us emphasize that exactly for these points in A4. the trivial fixed point 
g = g{x) = becomes unstable, for example, with respect to perturbations in g. 
To see this consider g = e: 

'^^ ae^ (17) 



rflnA 

where a < if (|T^ is true (i.e. arbitrary small e increases during the flow to low 
energies). This simple case illustrates the behavior of general system (p!2D. On the 
phase diagram it corresponds to // axis where the both RS and RSB phases exist. 

5 Summary 

Having started from the (space-time) supersymmetric Wess-Zumino model in a ran- 
dom and quenched background (|I|) we have found that renormalisation group equations 
(^ in a fixed point are quadratic homogenous equations in couplings and in g. The 
former property allowed us to take couplings very small as well as to reduce the moduli 
space to tM = RP"^. There are two types of points (phases) on this moduli space with 
either broken replica symmetry or not. 

Though we have found all IR-fixed points of one-loop /3-function, stability of nontriv- 
ial fixed points and analytic RG flow to them remain unexplored. Finally, it is interesting 
to generalize this analysis to more complex supersymmetric theories as well as to find 
realistic models whose critical behavior correspond to such theories. 
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